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We have studied the spin correlations with k= (|, |) propagation vector which appear be¬ 

low 0.4 K in Tb 2 Ti 207 spin liquid by combining powder neutron diffraction and specific heat on 
Tb 2 +ccTi 2 -cc 07 + 2 / samples with x=0, 0.01, -0.01. The k= (|, |, |) order clearly appears on all 
neutron patterns by subtracting a pattern at 1.2(1) K. Refining the subtracted patterns at 0.07K 
yields two possible spin structures, with spin-ice-like and monopole-like correlations respectively. 

Mesoscopic correlations involve Tb moments of 1 to 2 /xb ordered on a length scale of about 20 A. 

In addition, long range order involving a small spin component of 0.1 to 0.2 /xb is detected for the 
0 and 0.01 samples showing a peak in the specific heat. Comparison with previous single crystals 
data suggests that the (|, |, |) order settles in through nanometric spin textures with dominant 
spin ice character and correlated orientations, analogous to nanomagnetic twins. 


Rare earth pyrochlores with geometrical magnetic frus¬ 
tration are model compounds to investigate classical spin 
ice as well as quantum spin ice physics. The short range 
ordered ground states with large degeneracy give rise to 
a large variety of interesting phenomena such as mag¬ 
netic monopoleJI^^, first order transitiorP, or unconven¬ 
tional magnetoelastic excitationJ^. These frustrated 
ground states are highly sensitive to small perturba¬ 
tions and to the presence of defects. Already in the 
80 s, VillaiiP predicted that for an arbitrary small con¬ 
centration of defects the spinel lattice with pyrochlore 
structure could switch its ground state from spin liquid 
to spin glass. More recently many studies focus on the 
role of a minute amount of defects in pyrochlores. The 
crystal structure of R2Ti2 07 pyrochlores (where R^+ is 
a rare earth ion, Ti^+ is non magnetic) accommodates 
oxygen vacancies, small off-stoechiometries in the R^+or 
Ti^+ lattice, ’’stuffing” or site inversion. Oxygen vacan¬ 
cies introduce Ti^+ defects instead of Ti^+ ions, modify¬ 
ing the R^+ crystal field and inducing local distortions, 
with various consequences on the low temperature mag¬ 
netic properties and spin ice physics. In classical spin ice 
Dy2Ti207, oxygen depletion or a low level Dy stuffing 
in the Ti sites strongly slows down the monopole dy¬ 
namics at very low temperatur and affects the non¬ 
equilibrium macroscopic propertie d^ l ^^ l In the quantum 
spin ice Yb2Ti207, small off-stoechiometry clearly affects 
the magnetic ground stat^^. It appears that precise 
studies of the magnetic defects are needed to resolve con¬ 
troversies concerning these unconventional ground states. 

Among the pyrochlores, Tb2Ti207 is one of the most 
complex, the subject of numerous investigations since the 
discovery of its spin liquid character in 1999 . The Tb^+ 
short range correlated moments fluctuat^I^ down to 
20 mK at the time scale of the muon probe ( 10 “^s). The 


origin of these fluctuations has been highly debatecP^Hi^. 
Several mechanisms have been proposed taking into ac¬ 
count the non-Kramers nature of the Tb ion and the 
small energy gap between the ground and first excited 
crystal field doublets^which open the possibility of quan¬ 
tum spin ice behavioP^. Recent studies of the crystal and 
magnetic excitations using polarized neutrons show that 
the magneto-elastic coupling plays a crucial role, yield¬ 
ing an interaction between the first excited crystal field 
level and an acoustic phonon branch, resulting in a hy¬ 
brid phonon-magnon mod^^. The magnetoelastic cou¬ 
pling also induces original properties in th e pyrochlore 
family, such as giant magnet ostrictiorP^^ and very low 
thermal conduct ivit}dl. All these features suggest that 
the Tb2Ti207 ground state is determined by an interac¬ 
tion between quadrupolar moments, for which a dynam¬ 
ical Jahn-Tell er distor tion could be the first (mean-field) 
appr oximat iorP^^^^^^. 

In Tb2Ti207 an exotic transition has been observed 
around 0.4 K, as shown initially by anomalies of the spe¬ 
cific heafP^. Below this temperature, spin fluctuations 
coexist with s tatic spin correlations and spin glass-like 
irreversibilitie^^^^. In single crystals, these correlations 
have been detec ted by an energy analysis of the neu¬ 
tron cross sectiorP^^^^W^, yielding an elastic contribution 
with diffuse maxima at the positions of half integer Bragg 
peaks. This intensity vanishes under a small magnetic 
field of about 0.1 T and shows some sample dependence, 
suggesting the influence of a minute amount of defects. 
In powders, a systematic study of Tb2+ccTi2-a;07+y 
sample^^ showed that the system could be tuned by 
minute changes of the Tb content around a critical con¬ 
centration Xc=-0.005 from a purely spin liquid state for x 
< Xc to a state with partial antiferromagnetic order with 
k= (|, |) propagation vector for x > Xc- The detailed 
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nature of this k= (|, ^|) order was not determined in 
Ref. [3TJ but it seemed to extend over mesoscopic length 
scales, rather than the infinite length scale associated 
with symmetry breaking and Bragg scattering. 

To investigate this magnetic order, we studied pow¬ 
der samples of Tb 2 +a:Ti 2 -a: 07 +y with a precise control 
of the stoechiometry x= -0.01, 0 and 0.01, by combin¬ 
ing neutron diffraction and specific heat measurements. 
In all samples, by subtracting a pattern at 1.2(1) K, we 
observe magnetic peaks of Lorentzian line shape, which 
can be indexed with the k= (|, propagation vec¬ 

tor. These peaks correspond to the mesoscopic order 
discussed in the following. Its length scale deduced from 
the Lorentzian line shape is around 20 A, which corre¬ 
sponds to one magnetic cubic unit cell. Moreover, in the 
samples (x=0, 0.01) showing a transition in the specific 
heat, we also observe a Bragg-like contribution, involv¬ 
ing a minute Tb moment ordered at a length scale of 
about 500 A, close to the resolution limit. Both Bragg- 
like and Lorentzian contributions vanish around 0.4 K. 
By refining the magnetic patterns at the lowest temper¬ 
ature (0.07 K), we find two possible spin arrangements, 
fitting the data equally well. In both cases, the Tb mo¬ 
ments remain aligned along their local (111) axes, there¬ 
fore respect the Ising character of the crystal field. The 
first arrangement corresponds to magnetic moments in 
an ordered spin ice configuratiorP^, antiferromagnetically 
coupled along the (100) cubic axeJ^ without violating 
the ice rules. The second one is described by a piling 
along a [111] axis of monopole and antimonopole layers^ 
separated by spin ice layers. Comparing the calculated 
single crystal maps to the experimental one suggests that 
spin ice-like configurations with non-random mutual ori¬ 
entations, describing a ”nanomagnetic twin texture”, are 
actually favored at T = 0. The monopole-like struc¬ 
ture could be also stabilized when heating. We discuss 
the relationship between the mesoscopic order detected 
in the powder data and the diffuse scattering features 
(anisotropic excitations, pinch points) shown by single 
crystals. 


I. EXPERIMENT 

Powder samples of Th 2 -\-xT^h-x 07 -^y with x= 0, 0.01 
and -0.01 were synthesized as described in Ref. [311 The 
value of X was adjusted by changing the mass of the start¬ 
ing constituents, Tb 407 and Ti 02 , which were heated in 
air at 1350° C for several days with periodic grindings to 
ensure complete reactions. The obtained samples were 
ground into powder then annealed at 800° C in air for 
one day. Another stoechiometric sample, synthesized in 
2004 using a simila r procedurals and previously studied 
with neutron^^SMl g^go measured for comparison. 

All samples were characterized by X-ray diffraction 
and high resolution neutron diffraction patterns were also 
recorded for the stoechiometric samples on the spectrom¬ 
eter 3T2 at the Laboratoire Leon Brillouin (LLB), yield- 



FIG. 1: (color on line) Tb 2 Ti 2 07 . Neutron diffraction pat¬ 
terns measured at 0.07 K and 1.2 K, focusing on the region 
of the magnetic peaks, as shown by arrows. The full pattern 
at 0.07 K is shown in the inset, with the refinement of the 
crystal structure yielding the scaling of the ordered magnetic 
moment. 


ing a crystal structure in agreement with previous results. 
Slight colour changes between the samples with x=-0.01 
(beige), 0 (white), +0.01 (light brown) are attributed to 
small variations in the concentration of Ti^+ and oxygen 
ions induced by the off-stoechiometry. 

Powder neutron diffraction patterns were recorded on 
the diffractometer DIB at the Institut Lane Langevin 
(ILL) with a neutron wavel ength A = 2.52 A, down to 
0.07 K. As shown previousl}^^ a typical powder pat¬ 
tern of Tb 2 Ti 2 07 consists of nuclear Bragg peaks coex¬ 
isting with a large diffuse signal from the environmental 
background, incoherent scattering, and short range mag¬ 
netic correlations between first neighbors Tb^+ mome nts. 
These ” liquid-like” correlations start to appeaifl^l^ be¬ 
low about 50 K. The mesoscopic magnetic order studied 
here consists of small Lorentzian peaks which appear on 
the top of this diffuse signal on the pattern measured at 
0.07 K (Fig. [^. To isolate them from the other con¬ 
tributions, a neutron pattern measured at 1.2(1) K was 
subtracted. The amplitude of the ordered magnetic mo¬ 
ment was calibrated to that of the nuclear Bragg peaks 
by refining the crystal structure on the pattern at 0.07 K 
(Inset Fig. [^. 

Specific heat data Cp(T) were recorded below 20 K and 
down to 0.4 K using the relaxation method with a Phys¬ 
ical Property Measurement System (Quantum Design, 
Inc.). The phonon contribution to the total specific heat 
can be safely neglected below 8 K. In the temperature 
range of the measurements, i.e. above 0.4 K, the nuclear 
specific heat varies with the temperature as T~^. Follow¬ 
ing an earlier estimat^^, its contribution corresponds to 
0.35 and 0.155 J/K.mol at 0.4 and 0.6 K, respectively. 
The influence of the nuclear part in the upturn of the 
specific heat below 0.5 K for the x = +0.01 and x = 0 
samples is therefore marginal. 
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FIG. 2: (color on line) Tb 2 +ccTi 2 -cc 07 +^. (a) Total specific 
heat C'p(T) versus temperature T for several x contents. In 
inset the lattice constant of the cubic cell a deduced from X- 
ray data at 300 K (red dots), in comparison with values from 
Ref. [31] revised in 2015 (blue dots), (b) Cp(T) curves focusing 
on the low temperature region. In inset, the phase diagram 
of Ref. [31] where the arrows show the sample concentrations 
studied here. 
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II. RESULTS 

The lattice constant deduced from X-rays increases lin¬ 
early with X. The two stoechiometric samples have close 
values (a=10.15214(5) and 10.15229(5)1 at 300 K for the 
x=0 (2004) and x=0 sample respectively). The values for 
x=0 and 0.01 are slightly lower than reported in Ref. M 
by less than 5 10“^A (Inset Fig. |^). This difference 
could be connected with the initial characteristics of the 
starting constituents, such as the particle sizes of the 
starting powders, or with some differences in the ther¬ 
mal treatment procedures such as different cooling rates 
at the end of the thermal treatments and annealing. 

Specific heat data provide a direct evidence of the ex¬ 
otic transitiorP^. As shown in Fig. The Cp{T) curves 
of all samples overlap at high temperature. They strongly 
differ below 1 K, showing either a decrease (x=-0.01) or 
an upturn followed by a sharp peak (x=0, 0.01). These 
systematic changes with x reflect the phase diagram of 


FIG. 3: (color on line) Tb 2 +xTi 2 -cc 07 +y : Neutron magnetic 
patterns, focusing on the region of the most intense magnetic 
peaks. A pattern at 1.2(1) K was subtracted to extract the 
k= (1, 1, 1) mesoscopic order, (a) neutron patterns at 0.07 K 
for different Tb contents x. Intensities have been scaled to the 
intensity of a nuclear Bragg peak to allow a comparison be¬ 
tween the samples, (b) Neutron patterns versus temperature 
for the sample x=0 (2004). Short and long dashed lines shows 
the positions of (111) and (200) Bragg peaks respectively. 


Ref. |3T| The x=-0.01 sample situates in the ’’spin liquid” 
region of Ref. [311 For x=0 and x=+0.01, a phase transi¬ 
tion occurs at 0.4 K. The sample with nominal value x=0 
synthesized in 2004 does not show such an upturn, pos¬ 
sibly due to different characteristics of the starting con¬ 
stituents or slight differences in the thermal treatment. 
It confirms that a fine tuning of the off-stoechiometry re¬ 
quires a fully reproducible procedure in all the steps of 
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FIG. 4: (color on line) Tb 2 +ccTi 2 -cc 07 +y : Neutron magnetic 
patterns at 0.07 K for different Tb contents x. A pattern at 
1.2(1) K was subtracted. Scales have been enlarged and error 
bars omitted for clarity. Blue and red tick marks show the 
positions of the magnetic Bragg peaks for the mesoscopic and 
long range order respectively. Short and long dashed lines 
show narrow signals at the position of the (111) and (200) 
Bragg peaks respectively. Solid lines are the best refinements 
corresponding to spin ice-like and monopole-like structures, 
both yielding the same powder averaged intensity (see text). 

The magnetic patterns of all samples are shown in Fig. 
I^a) at the lowest temperature 0.07 K, focusing on the 
low angle region. Intensities were scaled to the intensity 
of a nuclear Bragg reflection to allow a direct compari¬ 
son between the samples. The patterns are quite similar, 
showing broad magnetic peaks with Lorenzian lineshape. 
In the (x=0, 0.01) samples, a small Bragg-like reflection 
of width limited by the experimental resolution can be 
detected on top of the strongest peaks. This Bragg-like 
contribution is clearly absent in the x=-0.01 sample. The 
Lorenzian peaks are indexed in the Space Group (SG) 
Fd3m with a k= (^, ^) propagation vector. The 

strongest | and (| ^ |) peaks are clearly seen at 
low angles. The (^ ^ |) peaks and weaker peaks with 
higher Miller indices can be detected at higher angles in 
the enlarged patterns (Figure 4). The temperature de¬ 
pendence of the magnetic pattern is shown in Figure 3(b) 
for the sample x=0 (2004). The mesoscopic magnetic 


order appears below about 0.3 K. Small narrow contri¬ 
butions are also detected at the positions of the ( 111 ) 
intense nuclear Bragg peak and at the (200) peak which 
corresponds to a nuclear extinction. There are respec¬ 
tively attributed to an imperfect data correction and to 
a minority phase with k=0 propagation vector similar 
to that of Tb 2 Sn 207 , and not considered in the analysis 
given below. 

III. NEUTRON DATA ANALYSIS 

We first performed a symmetry analysis in the SG 
Fd3m, with k= \) using the program Basireps 

of the Fullprof suit^^ and th e prog ram SaratP^. As al¬ 
ready discussed for Gd 2 Ti 2 0 '/^ 2 ®, the symmetry oper¬ 
ations which keep k invariant (symmetry operations of 
the small Group G/c) split the four rare earth sites of 

the pyrochlore lattice in two orbits: the site Tbi with 

local anisotropy axis along k remains invariant (orbit 1), 
whereas the three other sites Tb 2-4 are exchanged (orbit 
2 ). 

The decomposition of the magnetic representation Fxb 
for orbits 1 and 2 are: 

= ( 1 ) 

^ 2rF + (2) 

We show in Table [T| the Irreducible Representations 
(IR) and associated basis vectors (BV) from Sarah out¬ 
put, which provides real basis vectors. The four sites of 
a tetrahedron are Tbi (0,0,0) for the orbit 1 and Tb 2 
(|, |,0), Tbs (0, I) and Tb 4 (|,0, |), for the orbit 2, 
which are exchanged by the symmetry operations (xyz), 
(zxy) and (yzx). This table is similar to that given in 
Ref. ESI 

Systematic refinements using the Fullprof suite were 
performed to the data, for all possible combinations of 
the representations within the two orbits. A general 
state was constructed by combining the basis vectors for 
a given representation. Good refinements were obtained 
with F 3 and Tq for the orbit 1 and 2 respectively. Tak¬ 
ing into account the strong Ising anisotropy of the crys¬ 
tal field, we implemented additional constraints by con¬ 
straining the Tb moments to lie along their local <111> 
anisotropy axes and to have equal magnitude for the two 
orbits. With these constraints, the best refinement found 
from the symmetry analysis corresponds to the magnetic 
structure called ” monopole-like ” in the following. 

We also performed an extensive search in the cubic 
unit cell, of the magnetic structures compelling with the 
above constraints, searching for solutions starting from a 
Tb tetrahedron. Magnetic structures were refined using 
the Fullprof Suit^^. We found another magnetic struc¬ 
ture, called ” spin ice like” in the following, fitting the 
data equally well as the monopole-like structure. The 
corresponding refinements are shown in Figure As a 
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FIG. 5: (color on line) Tb 2 Ti 207 : Proposed magnetic structures for the Tb mesoscopic magnetic order. Arrows show the 
moment orientations, top: spin ice-like, bottom: monopole-like, a): primitive magnetic unit cells, b): cubic magnetic unit 
cells. Different colors correspond to different spin orientations, c): schematic magnetic stacking of the Tb tetrahedra along a 
[111] axis parallel to the propagation vector, as shown by a purple arrow. Blue and red tetrahedra correspond to monopole 
(one-in three-out) and antimonopole (one-out three-in) arrangements respectively. Green and orange tetrahedra correspond to 
spin ice (two-in two-out) arrangements, with magnetization respectively parallel and antiparallel to the b axis of the cubic cell. 
Dark grey tetrahedra correspond to spin ice arrangements with magnetization along a or c axes of the cubic cell. 


striking feature, these two spin structures have exactly 
the same powder averaged structure factor. They can 
be derived from each other by reversing the four mo¬ 
ments of appropriate Tb tetrahedra and performing an 
appropriate translation. Both correspond to complex non 
collinear antiferromagnetic structures, with cubic mag¬ 
netic unit cells of length 2a containing 128 Tb spins, in 
which the magnetization cancels. Drawings of these two 
spin structures are shown in Figure 

In the ’’spin ice like” magnetic structure (Figtop), 
the crystal cubic unit cell of length a contains 4 Tb 
tetrahedra with the same local “two-in two-out” spin ice 
structure. The magnetic cubic unit cell is obtained by 
translations along the (100) axes accompanied by spin 


reversal. This “spin-ice-like” structure, derived from the 
ordered spin ice structure of Tb 2 Sn 2 07 E 21 was previously 
proposed to account for the diffuse (^, ^) spin cor¬ 

relations in Tb 2 Ti 207 single crystaP^. It complies with 
the antiferromagnetic character of the magnetic unit cell 
by preserving the ice rules through ’’antiphase” walls in 
which the local magnetization of a Tb tetrahedron rotates 
by 90°. Relaxing the Ising condition from this structure 
yields a negligible canting angle (5 deg.) without notice¬ 
able improvement of the fit quality. We outline that the 
’’spin ice like” structure cannot be found by a symmetry 
analysis within the small Group G/c, since it violates the 
face centered cubic symmetry. It should be found within 
the decomposition of the General Representation in IR 
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TABLE I: Irreducible representations T and associated basis 
vectors ^ for the space group Fd3m with k= (|, |), cal¬ 

culated using the program SaralP^. The magnetic Tb atoms 
are at the sites Tbi: (0,0,0), Tb 2 : (|,|,0), Tbs: (0, |, |) 
and Tb4: (|,0, |). 


IR 

Basis Vector 

Atom 

BV 

components 




X 

y 

z 

Orbit 1 






Ta 

'01 

Tbi 

1 

1 

1 

Ts 

02 

Tbi 

1 

I 

0 


03 

Tbi 

1 

1 

2 

Orbit 2 







04 

Tb2 

1 

1 

0 



Tbs 

0 

1 

1 



Tb4 

1 

0 

1 


05 

Tb2 

0 

0 

1 



Tbs 

1 

0 

0 



Tb4 

0 

1 

0 

r4 

06 

Tb2 

1 

I 

0 



Tbs 

0 

1 

1 



Tb4 

1 

0 

1 

Te 

07 

Tb2 

1 

1 

0 



Tbs 

0 

I 

1 



Tb4 

1 

0 

1 


08 

Tb2 

0 

0 

1 



Tbs 

1 

0 

0 



Tb4 

0 

1 

0 


09 

Tb2 

0 

0 

0 



Tbs 

0 

1 

1 



Tb4 

1 

0 

1 


among the 4 branches of the k star, which could also 
provide multi-k structures. Such decomposition involves 
IR’s up to 12 dimensions, which is beyond the purpose of 
the paper. We also recall that multi-k structure cannot 
be distinguished from single-k ones in these powder data. 

The ” mo nop ole-like” structure (Fig bottom) found 
by symmetry analysis can be described in the same cu¬ 
bic unit cell as for the spin ice like structure by one 
monopole “one-in three-out” and three antimonopole 
“one-out three-in” tetrahedra. In other words, it con¬ 
sists of a piling along a [111] axis of monopole and anti¬ 
monopole layers in the Kagome planes, separated by spin 
ice layers. This “monopole-like” structure bears some re¬ 
semblance with the field induced antiferromagnetic one 
observed in Tb 2 Ti 207 under a high magnetic field ori¬ 
ented along a [110] axis, although the propagation vec¬ 
tor is differentP^. It complies with the antiferromagnetic 
character of the magnetic unit cell through a large num¬ 
ber of spin ice defects (namely static monopoles), pre- 


TABLE II: Tb 2 +xTi 2 -x 07 +j,: Mesoscopic magnetic moment 
M, domain size D and long range ordered magnetic moment 
Mlro. 


X content 

X — —0.01 

a: = 0(2004) 

X = 0 

X — +0.01 


0.86(1) 

2.17(3) 

1.89(1) 

0.65(1) 

d(A) 

21(2) 

24(2) 

22(2) 

22(2) 

T/lro(mb) 

0 

0 

0.2(1) 

0.15(3) 

C'p(T) peak 

No 

No 

Yes 

Yes 


serving the direction of the staggered magnetization. 

By using the program Isocif on the Bilbao Crystallo¬ 
graphic Serve]®!, we determine the magnetic space group 
and primitive unit cells (Fig[^) of the above structures, 
which are different from each other. The spin ice-like 
structure can be generated from a tetragonal magnetic 
unit cell of unit vectors um = = cl v^, cm = ‘^cl con¬ 

taining 64 Tb ions. The primitive magnetic cell of the 
monopole-like structure is monoclinic (/3 = 125.26°, um 
= a^=12.39A, bM = a^= 7.15A, cm = ay2=14.3lA) 
and contains 16 Tb ions. The magnetic space groups 
are of type IV, listed as follows: (BNS Ic-4c2 ; OG: P/- 
4’2m’; Litvin: 111.11.921) and (BNS: CcC ; OG: C 2 cm’; 
Litvin: 8.6. 43) fo r the spin ice and monopole structures 
respectivel 30 ®n. They are derived from the Fedorov 
space groups I-4c2 (NO. 120) and Cc (NO. 9). 

By fitting the neutron patterns to the above structures, 
we determined the ordered magnetic moment amplitude 
and the correlation length relative to the mesoscopic 
magnetic order for all samples (Table 2). These two 
quantities are the only free parameters of the magnetic 
patterns, and they can be determined separately. The 
amplitude of the ordered moment determines the peak 
intensities relative to the nuclear Bragg peaks, whereas 
the domain size controls the peak width. In all samples, 
the fitted correlation length of the (^, ^) correlations 

is around 22(2) A, close to the size of the cubic magnetic 
cell (20.3 A). The ordered moment is strongly enhanced 
for x=0, in the region of the critical point xc^ then de¬ 
creases both in the spin liquid and ordered regions. In 
addition, for the samples showing a peak of (7p(T), a 
minute long range ordered (LRO) moment of about 0.15 
/iB is detected. This LRO spin component orders over a 
length scale of about 500 A close to the resolution limit 
with the same magnetic structure. In all cases, the or¬ 
dered moments remain well below the local Tb mome nt 
of 5 /iB estimated from the crystal field schemd^^^ at 
low temperature. This means that in the ground state, 
a large component of the Tb moment is correlated only 
over first neighbor distances, both in the spin liquid and 
in the ordered region of the phase diagram. 

IV. DISCUSSION 

The above study focuses on mesoscopic (|, ^) corre¬ 

lations which emerge in Tb 2 Ti 207 powder patterns be- 
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FIG. 6: (color on line) Tb 2 Ti 207 single crystal: Elastic 

Diffuse scattering maps in a (hhO, 001) plane. The subtracted 
maps show the temperature evolution of the correlations, a): 
between 4K and 500 mK, the diffuse scattering starts to be 
structured, b): between 300 and 50 mK, the correlations do 
not evolve. 




low the exotic transition at T^', situated in the range 
0.3-0.4K, which presen ts some similarities with a spin 
glass transitiorP^H^^So]^ analysis provides a phe¬ 

nomenological quantitative description of these correla¬ 
tions, which involve only a fraction of the Tb moment, 
amounting to 0.1 to 2 //b depending on the length scale, 
namely well below the local Tb moment of 5 /tb- These 
correlations have been separated by subtracting a pattern 
above the transition. 

These findings should now be confronted to recent 
studies in single crystalJ^^^^^^. In single crystals, the 
temperature variation of the diffuse scattering shows a 
gradual evolution from a fluctuating spin liquid state 
characterized by isotropic first neighbor correlations at 
high temperatures, to a precursor regime (between T^ 
and 4K typically) where the diffuse scattering starts 
to be structured and the elastic component increases 
steadily, then to a regime (T< T^') where the elastic 
signal is dominant, well structured and T-independent. 


This temperature evolution is shown on the subtracted 
maps of Fig. We speculate that a similar evolution 
with temperature occurs in powders, so that the meso¬ 
scopic correlations can emerge. The precursor regime is 
not obviously seen here on the powder data due to the 
low signal background ratio, but it was actually seen in 
Tb 2 Sn 207 powdei!^, where the mesoscopic order occurs 
on a much larger length scale and involves a much higher 
fraction of the Tb moment. 

It has been previously arguecP^ that there are actually 
two transitions below 500 mK, the magnetic transition at 
Tg occurring about 0.15 K below the Cp{T) anomaly. 
This suggests two successive phenomena with decreasing 
temperature, related to quadrupolar and magnetic order 
respectively. In the ’’ordered” samples above Xc, the on¬ 
set of quadrupolar order seen by Cp{T) peak occurs at 
a higher temperature than the long range magnetic or¬ 
der, as frequently observed in rare earth inter met allicP^. 
In the ’’spin liquid” samples below Xc which show no 
anomaly of (7p(T), the magnetic order would remain lim¬ 
ited to short length scales. 
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FIG. 7: (color on line) Tb 2 Ti 2 07 : Diffuse scattering maps 
in a (hhO, 001) plane at 0.07 K. (a) : experimental map from 
Ref. m Galculated maps assuming a single magnetic domain 
of volume 8a^ with either spin-ice like (b) or monopole-like 
(c) structure; (d): Average over all single magnetic domains 
with equal probability. Monopole and spin ice -like structures 
yield the same average map. 


Below Tg the ela stic diffuse scattering measured on 
single crystalshows arm-like features along (111) 
directions and butterfly shaped features centered at (002) 
and (220). These features result in diffuse maxima 
around (^, ^) positions and equivalent (see Fig. [^). 

Moreover, the signal sharpens about certain wave vec¬ 
tors, defining pinch points in reciprocal space. In py- 


















rochlores, these pinch points typically result from the 
propagation of a local constraint govern ing the mutual 
orientations of the magnetic moments^^^^ such as the 
”two in/two out” rule in spin ices. The emergent gauge 
structure characterizes a Coulomb phas^^. Fennell et 
al have proposed that a similar ice rule should hold for 
Tb 2 Ti 207 , spin components along the local (111) axes 
being controlled by a ”two in/two out” condition, and in 
the transverse direction by a ”2-up/2-down” conditioiP^. 

To compare our phenomenological description with 
these features, we have calculated diffuse scattering maps 
assuming cubic magnetic domains of size 2 a and the spin 
orientations determined above, for different domain ori¬ 
entations (see the appendix for details). We compare 
these maps to those measured in the (hhO, 001) plane by 
elastic neutron scattering at 0.06(1) K shown in Figure 
First, we notice that the butterfly features around 
(002) or (220) as well as the pinch points which appear 
between the diffuse maxima in the experimental map of 
Fig[7^ are not reproduced by these calculations. Next, for 
the spin-ice-like structure, a few calculated maps, such 
as the map of Figure or the maps K1S2 and K1S3 in 
the Appendix, reproduce rather well the shape, extension 
and relative orientations of the diffuse maxima at (^, ^, 
|) and related positions. For the monopole-like struc¬ 
ture, the agreement is not as good, whatever the domain 
orientation (see an example in Figure [^). Averaging 
randomly over all domains orientations yields the same 
map for the two structures, shown in Figure [^, which 
is also in worse agreement with experiment. This sug¬ 
gests that specific spin-ice like textures, involving mag¬ 
netic domains of correlated orientations are preferred at 
the base temperature. They may be generated through 
nanomagnetic twins, induced by the lowering of symme¬ 
try from the Fd3m SG to the tetragonal or monoclinic 
SG’s discussed above. Such nano twins exist in crystal 
structureJ^ but their extension to magnetism is beyond 
the purpose of this study. 

To situate this picture in the context of a local con¬ 
straint, we notice that a local constraint supposes the 
existence of a manifold of degenerate configurations in 
the ground state. As temperature decreases, the sys¬ 
tem could select peculiar spin configurations out of this 
manifold. The spin-ice-like and monopole-like struc¬ 
tures determined above would belong to this manifold. 
This selection could be favored by small off-stoichiometry 
such as minute Tb substitution at the Ti sites that 
would slightly perturb the en ergy scheme by introducing 
further neighbor interact ions^^^^, local changes in the 
crystal fie ld, ra ndom magneto-elastic couplings or local 
distortionJ^^EH. These perturbations could reveal the ba¬ 
sic interactions by inducing correlated regions around im¬ 
purity sites, with the percolation of small long range or¬ 
dered moment above a certain level of substitution. Sim¬ 
ilarly, a low Fe or Go substitution of a Pd matrix reveals 
the high Pd band susceptibility by generating giant mo¬ 
ments around the Fe/Go impurity, with the onset of long 
range ferromagnetic order below one % substitutioiP^. 


One could wonder if such energetic selection supports 
low energy propagative excitations. This is the case in 
most fields of physics, where the stabilization of a meso¬ 
scopic order with static character allows exotic excita¬ 
tions to emerge. For instance, the slowing down of ther¬ 
mal relaxations allows dispersive excitations in super- 
paramagnetic particle^^, quantum tunneling in molecu¬ 
lar magnet^^, reversible magnetization in spin glasseJ^, 
or more trivially spin waves in classical magnets. Here we 
would expect excitations to stem from the k= (^, |) 

positions. Propagative excitations have been actually ob¬ 
served in Tb 2 Ti 207 but they stem from the pinch points 
at k=(lll). This qualitative difference suggests that 
these excitations are an intrinsic feature of the Coulomb 
phase. 

In summary, we have studied the mesoscopic order at 
play in Tb 2 +ccTi 2 -a^ 07 +y, which is tuned by small off- 
stoichiometry. We find the same type of correlations in 
all samples, with two possible structures, monopole-like 
and spin ice-like, describing the powder averaged mag¬ 
netic order. Magnetic moments of about 2 /tb are corre¬ 
lated over micro domains of 20 A size comparable to the 
cubic magnetic unit cell. Samples in the ordered region 
of the phase diagram involve in addition a very small 
long range ordered moment of less that 4% of the local 
Tb moment (^ 5 /ie)- Comparison with single crystal 
maps suggests that in the ground state, the spin-ice like 
structures are preferred and that they stabilized through 
twinned magnetic nano domains. Our observations pro¬ 
vide evidence that this mesoscopic order coexists in the 
ground state with a disordered Coulomb phase, which 
involves the major part of the local Tb moment. Their 
links remain however to be elucidated. 

We would like to acknowledge fruitful discussions and 
useful advices from J. Rodriguez-Carvajal, A. Wills, J. 
Robert, A. Gukasov, F. Damay, and A. Sazonov. We 
thank M. J. P. Gingras for a critical reading of the 
manuscript. We are grateful to the ILL for the use of neu¬ 
tron beam time and of the high-flux powder diffractome¬ 
ter DIB operated by the Centre National de la Recherche 
Scientifique (CNRS). 


V. APPENDIX: MAGNETIC DOMAINS 

In order to generate the magnetic domains for a given 
magnetic structure (either spin ice like or monopole like), 
we used the following procedure. For simplicity, we start 
from the cubic unit cell of lattice constant a containing 
16 Tb ions. The monopole and spin ice structure are 
defined within this cubic cell. We generate the K and 
S-domains within the SG Fd-3m. Among the 48 symme¬ 
try operations of the space group Fd-3m, the first 24 are 
associated with rotation matrices R with det(R)=l, and 
the 24 others deduced from them by an inversion center 
(det R=-l) do not need to be considered. Applying the 
24 symmetry operations to the propagation vector k= (^, 
|, |) splits them in 4 sub groups, which either leave k in- 
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variant or transform k in the k vectors 
|), (^, -^) of the star {k}. These four subgroups are 

related to the k-domains. For the little group Gk which 
leaves k invariant, the symmetry operations apply both 
to atomic positions and atomic moments, using the usual 
rules : i) classical spins transform as axial vectors under 
the symmetry operations; ii) for the symmetry operation 
which generate atomic positions out of the initial crys¬ 
tal cubic unit cell, equivalent atomic positions/magnetic 
moments in the cubic cell are obtained by appropriate 
translations Ri, such as ]\/Ejj — IN/Eqj exp ^2i7r k v^here 
Moj and Mij are the j moment of the initial cubic unit 
cell (1=0) and at the distance R^, respectively. This pro¬ 
cedure generates for each k-domain, 6 S-domains, each 
one corresponding to specific spin orientations of the Tb 
moments within the cubic unit cell. For each domain, 
the structure factor and neutron intensity are calculated 
in the (hhO, 001) plane, by taking the Fourier transform 
over a single magnetic cubic cell of size 2a, which corre¬ 


sponds to the typical correlation length deduced from the 
powder data. We then obtained 24 intensity maps, both 
for the spin ice-like and monopole-like structure. Figure 
shows these maps for the six S domains generated with 
G/c, where k= (^, ^). In the monopole structure none 

of the 24 maps corresponds to the single crystal data at 
low temperature, whereas in the spin ice structure some 
of them such as the K1S2 and K1S3 shown in Figure 
indeed bear some resemblance with the data. 

Performing an average over the intensities of the 24 do¬ 
mains with the same weight for each domain, yields the 
same intensity distribution for the spin ice and monopole 
structures (Figure [^) which does not well reproduce the 
experimental data. It suggests that specific spin tex¬ 
tures, associated with magnetic domains of correlated 
spin orientations, are actually stabilized. This ground 
state mesoscopic texture could consist of antiferromag- 
netically ordered spin ice blocks, with correlated orienta¬ 
tions of the staggered magnetization. 
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b) Monopole-like structure 
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FIG. 8: (color on line) Tb 2 Ti 2 07 : Diffuse scattering maps 49 
in the (hhO, 001) plane at 0.07 K, calculated for the six S- 
domains corresponding to the propagation vector k= (I? so 

1). Single magnetic domain of volume 8a^ with spin-ice like 
(a) and monopole-like (b) structures. si 
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